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Abstract 

The procedure of Dirac reduction of Poisson operators on submanifolds is discussed 
within a particularly useful special realization of the general Marsden-Ratiu reduc- 
tion procedure. The Dirac classification of constraints on 'first-class' constraints and 
'second-class' constraints is reexamined. 

AMS 2000 Subject Classification: 70H45, 53D17, 70G45 

1 Introduction 

Dirac bracket as well as Dirac's classification of constraints is nowadays a well recognized and 
very useful tool in the construction of Poisson dynamics on admissible submanifolds from a 
given Poisson dynamics on a given manifold. In this paper we consider the Dirac reduction 
procedure in a more general setting than is usually met in literature. In Section 2 we 
implement the Dirac reduction procedure into a particularly useful special realization of the 
general Marsden-Ratiu reduction scheme, based on the concept of transversal distributions. 
In Section 3 we reconsider the Dirac concept of first class constraints as it seems to be too 
restrictive. 

Firstly we recall few basic notions from Poisson geometry. Given a manifold Ai, a 
Poisson operator jt on M is a mapping n : T* Ai — > TA4 that is fibre-preserving (i.e. 
it\t*m '■ T*AA — > T x Ai for any x 6 Ai) and such that the induced bracket on the space 
C°°(Ai) of all smooth real- valued functions on Ai 

{.,.}„ :C 00 (M)xC 00 (M)^C QO (M) , {F, G}^ =' (dF, n dG) , (1) 
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where (.,.) is the dual map between TM. and T*M, is skew-symmetric and satisfies Ja- 
cobi identity (the bracket (JIJ always satisfies the Leibniz rule {F, GH}^ = G {F, H} n + 
H {F,G} n ). The symbol d denotes the operator of exterior differentiation. The operator 
7r can always be interpreted as a bivector, 7r £ A 2 (M) and in a given coordinate system 
(x 1 , . . . , x m ) on M we have 



^—f 9xj dxj 

i<3 ■> 



A function C : M — > K is called a Casimir function of the Poisson operator tt if for an 
arbitrary function F : M — > K we have {F, C}^ = or, equivalently, if irdC = 0. 



2 Marsden-Ratiu reduction for transversal distributions 

The Marsden-Ratiu reduction theorem £Q describes the procedure of reducing a Poisson 
operator n on arbitrary submanifold iS of our manifold M. This general procedure exists 
only if some conditions are satisfied. These conditions involve a distribution E (in the 
original notation of Marsden and Ratiu) that is a subbundle of TM. By a simple assumption, 
namely that this distribution is transversal, one can, however, satisfy all these conditions 
automatically. Below we reformulate the Marsden-Ratiu theorem in this more limited but 
useful setting. 

Consider an m-dimensional manifold M equipped with a Poisson operator tt and an 
s-dimensional submanifold S of M.. Fix a distribution Z of constant dimension k = m — s, 
that is a smooth collection of m-dimensional subspaces Z x C T X M. at every point x in M, 
which is transversal to S in the sense that no vector field Z S Z is at any point tangent to 
the submanifold S. Hence we have 

T X M = T X S © Z x 

for every x G S and, similarly, 

T*M=T*S®Z* X , 

where T*S is the annihilator of Z x and Z* is the annihilator of T X S. That means that if 
a is a one form in T*S then a(Z) — for all vectors Z g Z x and if [3 is a one-form in Z x 
then j3 vanishes on all vectors in TS X . 

Definition 1 A function F : A4 — > M is invariant with respect to Z if LzF = Z(F) = 
for any Z G Z. Similarly, a function F : M — > R is invariant with respect to Z on S 
(Z\g-invariant in short) if LzF\$ = Z(F)\g — for any Z G Z 

Here and in what follows the symbol Lz means the Lie derivative along the vector field 

Z. 

Definition 2 An operator n is called invariant with respect to the distribution Z if the 
functions that are Z- invariant form a Poisson subalgebra, that is, if F , G : M. —> K are 
two Z- invariant functions then {F, G}^ is again a Z- invariant function. Similarly, an 
operator 7r is called invariant with respect to the distribution Z on S if the functions that 
are Z\s - invariant form a Poisson subalgebra, that is, if F, G : M. — * M are two Z\$ - 
invariant functions then {F, G}^ is again a Z\s - invariant function. 



2 



We denote these Poisson subalgebras by A and As respectively. Let us observe, that if 
tt is iJl^-invariant for any manifold S u in a foliation of M. then it is also Z- invariant. 

Theorem 3 (Marsden and Ratiu [1]): LetS be a submanifold ofA4 equipped with a Poisson 
operator tt and let Z be a distribution in Ai that is transversal to S. If the operator tt is 
invariant with respect to the distribution Z on S, then the Poisson operator tt is reducible 
on S in the sense that on S there exists a (uniquely defined) Poisson operator ttr such that 
for any /, g : S — > R we have 

{f,9h R = {F,G}J s (2) 
for any Z\$ - invariant prolongations F and G of f and g respectively. 

The above construction, however, is difficult to perform in practice since it is often hard 
to find explicit expressions for the prolongations F and G. We now show how this difficulty 
can be omitted. 

Firstly, suppose that our submanifold S is given by k functionally independent equations 
tpi(x) = 0, i = 1, . . . , k (constraints) and that our transversal distribution Z is spanned by 
k vector fields Zi chosen such that the following orthogonality relation holds 

(d(fi, Zj) = Zj(ifii) = <%, (3) 

(this is no restriction since for any distribution Z transversal to S we can choose its basis 
so that is satisfied). We observe that in this case we have [Zi, Zj] tp^ = for all k, where 
[X, Y] — LxY — X(Y) — Y(X) is the Lie bracket (commutator) of the vector fields X, Y, 
so that [Zi,Zj] is always tangent to S. Then, in case that the distribution Z is involutivc 
(integrable), this means that [Zi, Zj] = for all Moreover, we define the vector fields 
Xi as 

Xi = ndtpi, i = 1, . . . , k. (4) 
There exists an important class of ^-invariant Poisson operators 

Lemma 4 If 

k 

L Zi TT = y £wf ) AZj i=l,...,k (5) 

3=1 

for some vector fields Wj*\ then the Poisson operator tt is invariant with respect to Z 
We sketch the proof here for the clarity of the text. 

Proof. Assume, that Lz t F = Lzfi = for all i. We have to show that L Zi {F, G} n = 
for all i, but, due to JSJ 

k 

L Zl {F, G} n = L Zi (dF, TrdG) = (dF, {wf A Z^dG) 

3=1 

since Lz^dF) — d(LziF) = (and similarly for G). On the other hand 

(dF, {Wf ] A Zj)dG^ = Z J {G)wf ) {F) - Z J {F)wf{G) = 
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since Zj{F) — L Zj F — (and similarly for G). ■ 

The condition is sufficient but not necessary. For example, if 



L Zt ir = ^yr j A[Z i ,Z j ] i = l,...,k 
i=i 

for some vector fields Wi, then the operator tt is also Z-invariant (one shows it by compu- 
tations similar to those in the above proof). In the case when tt satisfies JjJJ we apply the 
Lie derivative Lza to both sides of the equation Q. Due to JSJ we obtain 

[Z h Xi] = L Zj Xi = (Lzrfdtpi = w i U) A d Vi = 



(6) 



We observe that, if F and G are two ^-invariant functions and Vj are arbitrary 
vector fields, then (dF,J2jVj A Zj dG^ = since (dF, Vj A Zj dG) = Zj(G)Vj(F) - 
Zj(F)Vj(G) = on S. Thus the Poisson operator tt and its deformation of the form 

*D = tt - EjVj A Zj (7) 

in spite of the fact that they act differently on ^.5 both generate the same bracket on S so 
that both can be used to define our restricted operator 7Tr on S through @. Of course, 
the deformed operator ttd does not have to be Poisson, but nevertheless its restriction to 
S through (J2J) must be Poisson since it naturally coincides with similar restriction of tt to 
S. If we now consider a whole foliation of M. defined by the functions <fi with leaves S u (so 
that Sq = S) then it turns out that we can choose our (undetermined so far) vector fields 
Vj in |J7|) so that 

^D{a x ) £ T X S U for any a x G T*M and any x e M, (8) 
which has a far reaching consequence. 

Lemma 5 The deformation ttd given by that also satisfies |3J) is Poisson. 

Proof. The condition that ttd(c^x) is tangent to S u for any a x G T*M. is equivalent 
to the requirement that (dtfi, TTD(ct x )) = for all i. Due to the antisymmetry of ttd this 
requirement can be rewritten as (a x , irjj(ipi)) = for all i. Since a x is arbitrary, the 
condition attains the form irr)(dtpi) = for i = l,...,k. We now complete the set of 
functions <fi with some functions Xj to a coordinate system (2, ip) on M.. Then the matrix 
of the operator 7T£> has the last k rows and last k columns equal to zero while the m — k 
dimensional upper left block coincides on every leaf S u with the correponding ttr which is 
Poisson by the Marsden-Ratiu construction. ■ 
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Lemma 6 The condition can be written as 

k 

Vi-J2Vj(Vi) z j = x i- (9) 

3=1 

Proof. We know that the condition |JSJ can be written as nn(d(pi) — for i = 1, . . . , k. 
An easy calculation yields now that 

k 

= n D {dipi) = 7r(d(fi) - ^ ( z j(Pi)Vj ~ V j{fi) Z o) 

j 

k 

= x i -v i + Y,vAvi) z 3 

3=1 

due to the normalization condition ■ 

We now restrict ourselves to only two limit cases, when all Xi are tangent to S and when 
Xi span Z. 

2.1 The case when Xi are tangent to S 

We firstly assume that all the vectors X, are tangent to S. We have then naturally Xi(ipj) — 
0. This in turn means that {ipi,ipj} = (dipi,ndipj) = (dcpi,Xj) = so that all the vector 
fields Xi commute. In this case the simplest solution of © has the form Vi = X, and the 
corresponding deformation 10 attains the form 

k 

n D =TV-^X i AZ i . (10) 
i=i 

This deformation has been recently widely used for projecting Poisson pencils on symplectic 
leaves of one of their operators 0-0- 

Lemma 7 J3$ The vector fields in 0) can, in the case that all Xi are tangent to S , 

be chosen as tangent to S. 

Proof. Consider the projections Wj % of the vector fields wj onto S: 

r=l 

If are in Z, then = 0. The vector field Wj is indeed tangent to S since 

k 

W?\<p,) = Wf\ipi) - ]T wf{ Vr )5 lr = 0. 

r=l 
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Now 

k k k 

A Z, = £ ^Zj-J^ W f (<Pr)Z r A Zj 
3=1 3=1 3>=1 

the last term being equal to zero since L Zk {ipi,pj} n = implies Wj l \tp r ) = Wr\<p>j)- 

Thus Y* =1 Wf ] A Zj = £* =1 Wf A Zj. m 

Due to this gauge freedom, if we choose Wj as tangent to S (which means that 

wf\ip r ) = 0) then the formula ® yields that wf 1 = [Z h Xj] . Thus, due to the fact 
that we assumed 

k 

L Zt 7r = Y / [Z t ,X 3 ]AZ r (11) 

3=1 

Remark 8 In the case that the functions tpi are Casimir functions of n we have Xi — 
irdifi = so that the formula 111)) yields L Zi 7T — for all i, i.e. the vector fields Z% 
are symmetries of tt. In this case our reduction procedure 0) coincides with the standard 
projection onto a level set of Casimir functions (symplectic leaf in case there are no other 
Casimirs apart from (pi )Jffl/. 

From what we have said above it becomes clear that the above reduction scheme can 
be interpreted as a two-step procedure: firstly we deform the original Poisson tensor ir to a 
Poisson tensor ttd and then we obtain ttr as standard projection of ttd onto the level set S 
of its Casimirs (pi (thus we need not calculate the prolongations F and G in order to define 

Now we check what can be said about our vector fields Zi. 
According to Remark \^L Zi TTD = 0. On the other hand, due to IjlQI) . 

= L Zi tt d = J2 l Z i > Xj] A Zj - J2 L z, Xj A Zj — X, l\L z% Zj 

3 = 1 3 = 1 3=1 

so that Ylj=i Xj A [Zi, Zj] = 0. Of course one of the possible realizations of this condition is 
the case that the distribution Z be integrable since then [Zi, Zj] = 0. There are, however, 

other possibilities here. For example, if [Zi, Zj] = Yl s =i c tj-^ s w ^ n c tj = c sj' S J= i Xj A 
[Z h Zj] = as well. 

2.2 The case when Xj span Z 

This time we assume that Xi — ^2 k tpkiZk for some real valued functions (pij, which due to 
© yields 

<Pij = T,k<PkjZk(<Pi) = Xj(ipi) = {ifi,ipj} n . (12) 

The functions tpy define a fc-dimensional skew-symmetric matrix (p = {(pij) , i,j = 1, . . - k. 
The only condition imposed on ip is related to the demand that Xi span Z, i.e. det ip ^ 0. 
We thus do not have to assume JSJ this time since now the distribution Z is spanned by 
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the Hamiltonian vector fields Xi and thus ir is automatically invariant with respect to Z as 
Lxiit — for all i. It can be easily shown that 

[Xj,Xi] = X{ Vi , Vj }^ = 7T d{ip ll ip J } 7r = irdtpij. 

Now we look for solutions of 10 in the simple form Vi = aXi. Inserting this into 
and using the fact that tpy = — (fji we obtain 

k k 

= aX, - a^2Xj(tpi)Zj - X { = aX { + a^^ipjiZj - X l = (2a - l)X % 

3=1 3=1 

so that a = 1/2 and Vi = \Xi. In this case the deformation Q attains the form: 

k 

7r D =ir--^2x i AZ l (13) 

i=l 

and is, as mentioned above, Poisson. It is easy to check that our operator ttd defines the 
following bracket on M. 

k 

{F,G}^ D = {F, G}tt - {F^^A^hW^G}^, (14) 

where F, G : M. — > K are now two arbitrary functions on M, which is just the well known 
Dirac deformation |7] of the bracket {., .j^ associated with tt. 

Remark 9 // C : A4 —> R is a Casimir function of tt, then it is also a Casimir function of 
ttd since in this case {1$ yields 

m 

{F, C}t Td = {F, C}tt — ^2 {Fi ( Pi}Tr( ( P~ 1 )ij{ l Pj> C}tt =0 — = 0. (15) 

i.i=l 

We also know that the constraints (pi are Casimirs of the deformed operator ttd . Thus we 
can state that Dirac deformation preserves all the old Casimir functions and introduces new 
Casimirs (ft. 

It is now possible to restrict our Poisson operator tttj (or our Poisson bracket {., .},„) to 
a Poisson operator ttr (bracket {., -}it r ) on the submanifold S, i.e. the level set ipi = ... = 
ip m = of Casimirs of ttd , in a standard way. 

3 Existence of Dirac reduction 

We now present some realizations of the above Dirac case and discuss the classical concept 
of the Dirac classification of constraints. We will show that the classification of constraints 
as being either of first-class or of second-class, proposed by Dirac, should be reexamined 
when one looks at the problem from a more general point of view. 
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We recall that a constraint ipk is of first class if its Poisson bracket with all the remaining 
constants ipi vanishes on S, that is if 



0. 



(16) 



Otherwise tp k is of second-class. In the case that at least one of the constraints is of the 
first class, the matrix ipij in 112fl is singular on S so that the formula i|14fl cannot be used 
in order to define 7Tr. However, it may still be possible to define ttr via the above general 
scheme. This indicates that the concept of first class constraint is too narrow. Below we 
demonstrate the examples of Dirac reduction in case when constraints are of first class. 

We start with a simple example. Consider a 2n-dimensional manifold M. parametrized 
by coordinates (qi,...,qniPi,...,Pn) an d equipped with a Poisson operator of the form 





-Q, 



o 



where Q is a diagonal matrix of the form Q n — diag(qi i ... ! g n ). Consider a submanifold S 
given by a pair of constraints ipi(q,p) = q n = and f2(q, p) = Pn = 0. Then the matrix ip 
has the form 

q n 
-q n 



V = 



so that it is clearly singular on S (det(S') = on S) and 

q n 



o 

-l 



so that the Dirac formula (|14[) cannot be applied. However, the vector fields Z\ 
and Z2 = —q~ 1 Xi that span our distribution Z are not singular on S since X\ = - 
and X2 — q n d/dq n so that the deformation 1)1 3[) becomes 



= qn'x 2 

q n d/dp n 



1 v v d 

ir D = it Xi A X 2 = 7r - q n - — A - — 

q n aq n dpn 



\ - a 

dqi dpi 



i=i 



and is clearly reducible on iS. The operator itr obtained on S parametrized by coordinates 
(qi,...,q n -i,Pi,...,Pn-i) is 

Q„_i 

-Qn-l 

This simple example clearly illustrates that Dirac's classification is too strong. As a second 
example we consider a particle moving in a Ricmannian manifold Q of dimension three with 
a contravariant metric tensor 

"001 

G= 1 

1 

given in some coordinates (g 1 , <? 2 , g 3 ). Suppose that this particle is subordinated to a holo- 
nomic constraint on Q given by 



<pi(q) = gY + g 3 = 0. 



(17) 
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This defines a submanifold of Q. The velocity v — Xw=i v l d/dq l of this particle must then 
remain tangent to this submanifold so that 

i—l 

and thus in our coordinates v l — Y^j G^pj the motion of the particle in the phase space 
M. = T*Q is constrained not only by l|17|l but also by the relation 

(p 2 (q,p) = J2 G ij ^f Pj =Pl +P2q 1 +P3Q 2 = (18) 

that is nothing else than the lift of (|17|l to M.. The constraints (|17fl - l|18|) define a four- 
dimensional submanifold S of M.. We now introduce the following Poisson structure on 
M: 












q 1 


-1 
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q 





-1 











2q 3 


q 2 


q 1 


-q 1 


-q 2 


-2q 3 





P2 


Pa 


1 





-q 2 


~P2 











1 


-q 1 


-P3 









Again the matrix if is singular, since c/?i2 = 2(q 1 q 2 + q ) = 2(pi which obviously vanishes 
on S. One can, however, perform the deformation (|13|l . A quite lengthy but straightforward 
computation shows that in this case 



tt d = 












q 1 


-1 














q 2 





-1 











-2q x q 2 


q 2 


q 1 


-q 1 


-q 2 


2q 1 q 2 





P2 


P3 


1 





-q 2 


-P2 











1 


-q 1 


-P3 









and this operator can be projected onto S. To do this, one can first pass to the Casimir 
variables 

(q 1 , q 2 , (q) , (q, p),p%,P3) 

since, due to the fact that it is easiest to eliminate q 3 and p\ from the system of equations 
(pi = (pi(q) = 0, f2 = (P2{q,p) = 0, we parametrize our submanifold by the coordinates 
(g 1 , q 2 ,P2,P3)- In these variables the operator 7Tr attains the canonical form 









-1 














-1 


1 














1 









Our two examples show that the condition 116fl is only a necessary condition for nonex- 
istence of 7Tr on iS, but is not a sufficient one. Hence the definition of first class constraints 
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has to be made weaker. Even in the case when we deal with a real first class constraint we 
can obtain ttr on S coming from the Dirac reduction of n. We demonstrate this below. 

Firstly we assume that we have a pair of constraints (pi , ip2 that define our submanifold 
S = {(pi = 0, ip2 = 0} and such that they are of second class, i.e. that ^12 |s = {<Pi, ^2} \s 7^ 
0. It is clear that our submanifold S can be parametrized in infinitely many different ways 
by constraints tpi = 0, tp 2 = 0, where 

(pi = IpHPl + V>2V2, (f>2 = VWl + "04^2 (19) 

and where ipi are some functions on M. such that ipi\s 7^ and such that 

D((pi,(p2) 



D = 



D{(pi,(p 2 ) 



^1^4 - Ms + 0. (20) 



One can prove the following 



Lemma 10 The deformations $13p given by the pair <pi,v?2 of constraints and by the pair 
(pi, (f>2 of constraints define the same reduced Poisson operator ttr on S. 

Proof. For the moment we denote the deformation l|13fl defined through (pi,(p 2 by ttd 
and the corresponding deformation defined through (fi, <p2 by ttd. Applying (|13fl we easily 
get that for any two functions A, B : M. — > R 

\A B\ =\A B\ + i- 8 ' ~ {4 Pihr t 8 . 

where we have assumed that {(pi, ^2}^ does not vanish on S. Similarly 

{a, B ho = { A, By + {AMAM^^iMAM^k , (21) 

where {<pi,<p2}„ does not vanish on S due to 12U|) . Using the relations (|19l) between the 
deformed constraints (fn and the original constraints (pi, the Leibniz property of Poisson 
brackets and the fact that the functions (pi vanish on S we obtain 

{(PI,(P2}J S = D {<PU<P2}J S 

and 

({A, {B, $i}„ - {A, (pi}„ {B, fc} v )\ s 

= D ({A, V 2\« {B, <pi} v - {A, ^ {B, (p 2 }J\ s 

so that the nonzero terms D in the numerator and denominator of l|21|) cancel and we obtain 
{^4, B} 7ro I = {A, B}~ d I which implies that the projections of ttd and ttd onto S coincide. 
■ 

In this nonsingular case the distribution Z along which we project a Poisson tensor ir 
usually changes after reparametrization, but Z\$ remains the same as can be easily demon- 
strated. Thus in case of the second class constraints one has a "canonical" way of projecting 
7r onto S. 
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We now suppose that the constraints ifi are of first class, that is {<pi, tp2}n\s = and 
that the singularity in ttd is not removable. We may still attempt to define the projection ttr 
by reparametrizing S as in l|19fl above. It turns out that among an infinite set of admissible 
reparametrizations there are some exceptional which, although they fulfil the condition 
jnjj, nevertheless eliminate the singularity in ttd- In this case, however, by choosing a new 
parametrization tf>2 of S we change the distribution Z even on S so that we cannot 
expect that the projection ttr will be independent of the choice of the parametrization. We 
lose a natural, "canonical" choice of projection, but we still can perform the projection, 
although in infinitely many nonequivalcnt ways. We illustrate this below in a sequence of 
examples. 

Consider a six-dimensional manifold M. parametrized with coordinates (91,92, 93jPi;P2,P3) 
with the following Poisson operator: 



7T = 












1 


9i 














9i 


292 + 1 


93 














93 





-1 


-91 








-pl 





-91 


-292 - 1 


-93 


Pl 





P3 





-93 








~P3 






Consider now a four-dimensional submanifold S in Ai given by the relations 

Vi(<l,P) = © = 0, (f 2 (q,p) =P3 = 0. (22) 

It is clear that {ipi, ^2}^ vanishes on the whole manifold Ai (and thus on S) so that these 
constraints do not define any Dirac deformation at all. We now deform \\Tll as 

<Pl = fl + <P2, <?2 = (~P2 ~ qiPl)fl + <f2 (23) 

Calculation shows {ipi, ^2}^ — {P3 — 93)93 so that {^1, <^2}w|s = 0. One can show that after 
introducing the Casimir variables (qi, q%, <pi,Pi,p%, $2) the deformed operator ttd attains 
the form 








9193 

93 -P3 





1 


-91 


" 




O 91 13 





91 - 


f 2 Piqa 


-ql + 9 







13 -P3 


93 ~P3 
























-1 


9! 2 Pl93 

* 93 —P3 








-pi 







9i 


- q l-e 





Pl 





























where now 93 = 93(9, p, tp) and p 3 = p 3 {q,P, f) and 9 = (93 + p 2 q 3 + 9i93Pi)/(93 - P3), and 
as such is clearly singular on S and thus unreducible. This situation seems to be the most 
common, i.e. a spontaneous choice of parametrization almost always leads to a singularity. 
However, if we perform a slightly different deformation of (1221 : 

fit = tpi, <^2 = {~P2 - qiPi^Pi + V2 (24) 
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so that {<pi, V?2jv = —% is again zero on S, then the operator irp becomes nonsingular and 
its projection on S has the following form 






-1 -9i 



9i 




-9i 
-if 
Pi 




9i 9i - 1 -Pi 

in the variables (91, 92, Pi, P2)- Yet another deformation (even this time of the form (|19|l '): 

¥>1 = 92<y5l , <f2 = (P2 + </?2 Vl (25) 

yields a quite complicated expression on {(pi, <^2jv 

{PuPs}* = (392 + 1)93 + 93, 

so that it again vanishes on S, but ttd is again nonsingular and in the same variables 
(9i? 92, Pi; P2) its projection becomes 



3<;2 + l 





- 1 






91 92 
392 + 1 





1 



3q 2 + l 
9192 



3(J2 + 1 



91 P2 

392+1 






91P2 

3<J2 + 1 





which concludes our series of examples. 



4 Conclusions 

In this article we have focused on two issues involving Dirac reductions of Poisson operators 
on submanifolds. In the first part of the article we have shown how the Dirac reduction 
procedure fits in a natural way, i.e. as a result of two natural assumptions about the defor- 
mation 7Td of 7r, in the general Marsden-Ratiu reduction scheme. In the second part of our 
considerations we have demonstrated that the Dirac reduction procedure is often possible 
even in cases when the constraints that define our submanifold are of first class (in Dirac 
terminology), possibly after some suitably chosen reparametrization of the submanifold S. 

Acknowledgement One of the authors (M.B.) would like to thank Department of Science 
and Technology, Linkoping University, for its kind hospitality during his stay in Sweden in 
the fall of 2002. 
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